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Abstract 



Dark matter caustics have specific density profiles and, therefore, precisely 
calculable gravitational lensing properties. We present a formalism which 
simplifies the relevant calculations, and apply it to four specific cases. In the 
first three, the line of sight is tangent to a smooth caustic surface. The cur- 
vature of the surface at the tangent point is positive, negative or zero. In the 
fourth case the line of sight passes near a cusp. For each we derive the map 
between the image and source planes. In some cases, a point source has mul- 
tiple images and experiences infinite magnification when the images merge. 
Unfortunately, for the dark matter caustics expected in realistic galactic halo 
models, the angular resolution required to resolve the multiple images is not 
presently achievable. A more promising approach aims to observe the distor- 
tions caused by dark matter caustics in the images of extended sources such 
as radio jets. 
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I. INTRODUCTION 



Gravitational lensing techniques have proven useful in studying the distribution of dark 
matter in the universe. They have been used to reveal the existence of massive compact 
halo objects (MACHOs) in galaxies [1], and to constrain the mass distribution in galaxy 
clusters [2]. In this paper, we calculate the gravitational lensing properties of dark matter 
caustics [3-5]. Gravitational lensing by dark matter caustics has been discussed by Hogan 
[6]. We confirm Hogan's results for the case he considered, called the 'concave fold' in 
our nomenclature. We study additional cases and introduce a formalism to facilitate the 
calculations. 

In the cold dark matter (CDM) cosmology, caustics in the distribution of dark matter 
are plentiful once density perturbations enter the non-linear regime. A dark matter caustic 
is a surface in physical space where the dark matter density is large because the sheet on 
which the particles lie in phase-space has a fold there. Dark matter caustics can move about 
but they are stable. 

To see the validity of the preceding statements, note that before the onset of galaxy 
formation the cold dark matter particles lie on a thin 3-dim. sheet in 6-dim. phase-space. 
The thickness of this sheet is the primordial velocity dispersion Sv of the particles. It is of 
order [4] 



la ) 



M^)-3-10-M^^) [^] (1.1) 



for axions, and 



. 10- (^)"' (^f (1.2) 



for weakly interacting massive particles (WlMPs), such as the lightest supersymmetric par- 
ticle in supersymmetric extensions of the Standard Model. In Eqs. (1.1-2), is the present 
age of the universe and nia and mw are respectively the masses of the axion and WIMP. 
The present average axion number density is 



1.5 10^4 „ ( hV /lO'^eV^ 



= _ Si„ (1.3) 

where Q,a is their present energy density in units of the critical density, and h is the Hubble 
expansion rate in units of 100 km/(s-Mpc). Likewise, for WIMPs 



1.5 10^" „ h 
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pc^ \mw J 

The large exponents in Eqs. (1.3-4) indicate that the density of cold dark matter particles is 
enormous in terms of astronomical length scales. Being effectively collisionless, the particles 
experience only gravitational forces. These are universal and vary only on huge distances 
compared to the interparticle distance. Hence the sheet on which the particles lie in phase- 
space is continuous. It cannot break and therefore its evolution is constrained by topology. 
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In contrast, in present numerical simulations of galactic halo formation, the particle mass 
m is of order 10^ Mq and hence the number density 



The resolution of the simulations is inadequate to describe the dynamical evolution of the 
phase-space sheet and the physical implications of its existence. 

Where a galaxy forms, the sheet wraps up in phase-space, turning clockwise in any two 
dim. cut {x, x) of that space, x is the physical space coordinate in an arbitrary direction 
and x its associated velocity. The outcome of this process is a discrete set of flows at any 
physical point in a galactic halo [7] . Two flows are associated with particles falling through 
the galaxy for the flrst time (n = 1), two other flows are associated with particles falling 
through the galaxy for the second time (n = 2), and so on. Scattering in the gravitational 
wells of inhomogeneities in the galaxy (e.g. molecular clouds and globular clusters) are 
ineffective in thermalizing the flows with low values of n. 

Caustics appear wherever the projection of the phase-space sheet onto physical space has 
a fold [3-5]. Gcncrically, caustics are surfaces in physical space. They separate regions with 
differing number of flows. On one side of a caustic surface there arc two more flows than on 
the other. Because caustic surfaces occur wherever the number of ffows changes, they are 
topologically stable, in direct analogy with domain walls. At the caustic, the dark matter 
density is very large. 

There are two types of caustics in the halos of galaxies, inner and outer. The outer 
caustics arc simple fold {A2) catastrophes located on topological spheres surrounding the 
galaxy. They occur where a given outffow reaches its furthest distance from the galactic 
center before falling back in. The inner caustics are rings [3]. They are located near where 
the particles with the most angular momentum in a given inflow reach their distance of 
closest approach to the galactic center before going back out. A caustic ring is a closed tube 
whose cross-section is a D_4 (also called elliptic umbilic) catastrophe [4]. The existence 
of these caustics and their topological properties are independent of any assumptions of 
symmetry. 

Dark matter caustics have very well deflned density proflles, and hence calculable grav- 
itational lensing signatures [6]. It is the purpose of this paper to derive these signatures 
in a number of speciflc cases. In the limit of zero velocity dispersion {Sv = 0), the density 
diverges when one approaches a caustic surface, on the side which has two extra flows, as 
the inverse square root of the distance to the surface. This divergence is cut off if there 
is velocity dispersion, because the location of the caustic surface gets smeared over some 
distance 5x. For the dark matter caustics in galactic halos, Sx and Sv are related by [3] 



where v is the order of magnitude of the velocity of the particles in the flow and R is the 
distance scale over which that flow turns around, i.e. changes its direction. For a galaxy 
like our own, v = 500 km/s and R = 200 kpc are typical orders of magnitude. 

As was mentioned earlier, the primordial velocity dispersion of the leading cold dark 
matter candidates is very small. Using Eqs. (1.1-2), one flnds that axion caustics in galactic 
halos are typically smeared over 
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as a result of their primordial velocity dispersion, whereas WIMP caustics arc smeared over 



It should be kept in mind, however, that a cold dark matter flow may have an effective 
velocity dispersion which is larger than its primordial velocity dispersion. Effective velocity 
dispersion occurs when the sheet on which the dark matter particles lie in phase-space is 
wrapped up on scales which are small compared to the galaxy as a whole. It is associated with 
the dumpiness of the dark matter faUing onto the galaxy. The effective velocity dispersion 
of a flow may vary from point to point, taking larger values where more small scale structure 
has formed, and taking the minimum primordial value where no small scale structure has 
formed. For a coarse-grained observer, the dark matter caustic is smeared over 5x given by 
Eq. (1.6) where 5v is the effective velocity dispersion of the flow. 

Little is known about the size of the effective velocity dispersion of dark matter flows 
in galactic halos. However, in ref. [8], a triangular feature in the IRAS map of the Milky 
Way was interpreted as the imprint upon baryonic matter of the caustic ring of dark matter 
nearest to us. The sharpness of the feature's edges implies an upper limit of 20 pc on the 
distance 5x over which that caustic is smeared, and hence an upper limit of order 50 m/s 
on the effective velocity dispersion of the corresponding flow. 

The gravitational lensing effects of a caustic surface are largest when the hue of sight 
is near tangent to the surface because the contrast in column density is largest there. The 
effects depend on the curvature of the caustic surface at the tangent point in the direction 
of the line of sight: the smaller the curvature, the larger the effects. A caustic is an oriented 
surface because one side has two more flows than the other. We will consider three cases of 
gravitational lensing by a smooth caustic surface. In the flrst case, the line of sight is near 
tangent to a caustic surface which curves towards the side with two extra flows; see Fig. 4. 
We call such a surface 'concave'. In the second case, the surface is 'convex', i.e. it curves 
away from the side with two extra flows; see Fig. 6. In the third case, the caustic surface 
has zero curvature at the tangent point (the radius of curvature is inflnite), but the tangent 
line is entirely outside the side with two extra flows. 

Caustic surfaces may have cusps. The outer dark matter caustics of galactic halos are 
topological spheres, which have no cusps. But the inner dark matter caustics of galactic 
halos are closed tubes whose cross-section is a D_4 catastrophe; see Fig. 2. This cross- 
section has three cusps. We call it a 'tricusp'. The fourth case of gravitational lensing which 
we consider has line of sight near a cusp, and parallel to the plane of the cusp; see Fig. 9. 

Gravitational lensing produces a map of an object surface onto an image surface. The 
magnification is the inverse of the Jacobian of this map. Because dark matter caustics have 
well defined density profiles, it is a neat mathematical exercise to calculate their gravitational 
lensing characteristics. The images of extended sources may show distortions that can 
be unambiguously attributed to lensing by dark matter caustics in the limit of perfect 
observations. We will see that in three of the cases considered, a point source can have 
multiple images. In those cases, when two images merge, the Jacobian of the map vanishes 




(1.8) 
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and the magnification diverges. So, at least in theory, it seems that gravitational lensing is 
a very attractive tool for investigating dark matter caustics. Observation of the calculated 
lensing signatures would give direct evidence for caustics and CDM. 

We have been particularly motivated by the possibility [6] that the observer might be 
able to distinguish between axions and WIMPs by determining the distance over which the 
caustics are smeared. The nearby caustic, whose position is revealed according to ref. [8] by 
a triangular feature in the IRAS map of the Milky Way plane, is only 1 kpc away from us in 
the direction of observation. By observing the gravitational lensing due to that caustic, one 
may be able to measure 5x as small as 10^^ cm, assuming an angular resolution of 3 • 10~^ 
radians. If 6x turned out to be that small, the WIMP dark matter hypothesis would be 
severely challenged; see Eq. (1.8). Unfortunately, as will be shown below, the gravitational 
lensing due to a caustic only a kpc away from us is too weak to be observed with present 
instruments. It is well known that gravitational lensing effects are proportional to ^^^^^ 
where Ds^ Dl and Dls are respectively the distances from the observer to the source, from 
the observer to the lens and from the lens to the source. We will see below that, for the 
gravitational lensing effects of dark matter caustics to be observable with present technology, 
the lenses and sources must be as far as possible, at the cosmological distances of order Gpc. 
Even then, the observation of such effects will be difficult. Unfortunately, at Gpc distances 
it is not possible to measure 5x as small as Eqs. (1.7-8) with foreseeable technology. So 
it seems unlikely that one will be able to distinguish between dark matter candidates on 
the basis of the gravitational lensing characteristics of the caustics they form. Henceforth, 
unless otherwise stated, the velocity dispersion is set equal to zero. 

The remainder of this paper is organized as follows. In section II, we describe the outer 
and inner caustics of dark matter in galactic halos. We provide estimates of their sizes 
and of the strengths of their density profiles, using the self-similar infall model of galactic 
halo formation. In section III, the general formalism of gravitational lensing is reviewed. 
We show how the calculations can be streamlined for the case of gravitational lensing by 
dark matter caustics. In section IV, we calculate the gravitational lensing properties of 
dark matter caustics in the four cases mentioned above. In section V, we summarize our 
conclusions. 

II. DENSITY PROFILES OF CAUSTICS 

In this section we discuss the density profiles of the dark matter caustics present in 
galactic halos. We consider the outer caustics first. 

A. Outer caustics 

Outer caustics are closed surfaces (topological spheres) near the (n + l)th turn-around 
radii with n — 1, 2, 3, ... . Indeed the number of flows changes by two there because of the 
fall back of the particles. There is no caustic associated with the first turn-around, see Fig. 
1. The outer caustics are described by simple "fold" (^42) catastrophes. Their density profile 
is: 
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d» = ^e(a) (2.1) 

for small cr, where a is the distance to the caustic, 0(cr) is the Heaviside function, and A„ 
is a constant which we call the jold coefficient, o" > on the side with two extra flows, i.e. 
towards the galactic center. Therefore, when an outer caustic is approached from the inside 
the density diverges as ~ abruptly falling to zero on the outside. The observation [9] of 
arc-like shells surrounding giant elliptical galaxies can be interpreted [10] as outer caustics 
in the distribution of baryonic matter falling onto those galaxies. 

To estimate An in Eq. (2.1), consider the time evolution of CDM particles which are 
falling out of a galactic halo for the nth time and then fall back in. Let Rn be the turn- 
around radius. We assume that the rotation curve of the galaxy is flat near r = i?„ with 
time-independent rotation velocity Vrot- The gravitational potential is then: 

nr)=^rltln(-^) . (2.2) 

The particles have trajectory r(t) such that 

dr 



dt 




2{E-V{r))^v,,,\2\n{^] , (2.3) 



where E is the energy per unit mass. Eq. (2.3) neglects the angular momentum of the 
particles. This is a good approximation at turnaround since the particles are far from their 
distance of closest approach to the galactic center. Finally, we assume that the flow is 
stationary. In that case the number of particles flowing per unit solid angle and per unit 
time, is independent of t and r, and the caustic is located exactly at the (n -|- l)th 

turn- around radius 

Let us emphasize that none of the assumptions - spherical symmetry, flat rotation curve, 
time independence of the gravitational potential, radial orbits, and stationarity of the flow - 
affect the existence of the outer caustics or the fact that they have the density profile given 
by Eq. (2.1). The assumptions are made only to obtain estimates of the coefficients An. 

The mass density of particles, d,n{r), follows from the equality: "^-^^dt = ^^^^r'^dr\ 
where m is the mass of each particle. The factor of 2 appears because there are two distinct 
flows, out and in. Using Eq. (2.3) we obtain the density distribution near the nth outer 
caustic: 



dn{r) = 2- 



dndt 



n rVotJ21n(f^) 



(2.4) 



where = rn-^^. Near the caustic. In (^^^ = where a — Rn — r. Inserting this into 
Eq. (2.4) and comparing with Eq. (2.1), we find: 



V2 dM 

A, 



Vrot dVtdt 



K'" ■ (2.5) 



n 



dM 



Estimates of i?„ and ^^^^ 
[12,13]. The infall is ca 



can be extracted from ref. [11] for the case of self-similar infall 

n 

led self-similar if it is time-independent after all distances are 
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rescaled by the turn-around radius R{t) at time t and all masses are rescaled by the mass 
M{t) interior to R(t). In the case of zero angular momentum and spherical symmetry, the 
infall is self-similar if the initial overdensity profile has the form = ("m^) where Mq and 



e are parameters [12]. In CDM theories of large scale structure formation, e is expected to 
be in the range 0.2 to 0.35 [11]. In that range, the galactic rotation curves predicted by the 
self-similar infall model are flat [12]. 

The Rn and do not depend sharply upon e. Our estimates are for e = 0.2 because 

they can be most readily obtained from ref. [11] in that case. For e = 0.2, the radii of the 
(n + l)th turn-around spheres are approximately 

{Rn:n = 1, 2, ...} ~ (240, 120, 90, 70, 60, ...) kpc • ( • (2.6) 



Moreover, one has 



dM 



dCtdt 



Z^ = ^nx/2^ , (2.7) 



with 



{F„ : n = 1, 2, ...} ~ (20, 8, 5, 4, 3, ...) 10"^ . (2.8) 
Combining Eqs. (2.5 - 2.8), we find 



{A„:n=l,2,..}~(2, 2, 2, 3, 3,...) 



10 ^gr 
cm^ kpc^/^ 

\ 1/2 / , \ 3/2 



220km/sy VO.7^ 



(2.9) 



Generally the outer caustics are 'concave', i.e. they are curved towards the side which has 
two extra fiows. Their radius of curvature is of order R^. The lensing by concave caustic 
surfaces is discussed in section IV.A. 

In the next subsection, we discuss the density profiles of inner caustics. 



B. Inner caustics 

The inner caustics [3] are closed tubes whose cross-section, shown in Fig. 2, is a D_4 
catastrophe [4] . They are located near where the particles with the most angular momentum 
in a given infiow are at their distance of closest approach to the galactic center. We call them 
'caustic rings'. For simplicity, we study caustic rings which are axially symmetric about the 
z-direction as well as refiection symmetric with respect to the z — plane. The dark matter 
flow is then effectively 2-dimensional. In galactocentric cylindrical coordinates, the flow at 
such a caustic ring is described by [4]: 

z{a,T) = bar and p{c(,t) = a + —{r — tq)'^ — -a^ , (2.10) 
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to lowest order in an expansion in powers of a and r. a and r are parameters labeling the 
particles in the flow, p is distance to the 2;-axis. r is the time a particle crosses the z = 
plane, a is the declination of the particle, relative to the z = plane, when it was at last 
turnaround, b, a, u, tq and s are constants characteristic of the caustic ring, b, u, s and 
a are positive, whereas tq can have either sign, a is the radius of the ring; see Fig. 2. The 
transverse dimensions of the ring are p = ^uTq and q — in the p and z directions 

respectively. The cross-section of the ring has three cusps. So we call it the "tricusp". 
The cusps are topological properties of the caustic. Their presence does not depend on any 
assumptions of symmetry. 

The physical space density is given by [4] 



, , , 1 " dM . , cos (a) 



p . dfldr 



|^2(«,r) 



(2.11) 



where 



dM 



dM 



^^^^ is the mass falling in per unit time and unit solid angle, aj (p, z) 
l,...,n, are the solutions of p{a,T) — p and z{a,T) — z. n{p,z) 



dQ,dT 2iT cos (a) 

and Tj{p,z), with 

is the number of flows at {p,z). n — 2 outside the caustic tube, whereas n = 4 inside. 
D2{ot,T) = det D(Q;,r) with 



I d2_ 
' da 



dz 
da 



1 9p d£ 
^ dr dr 



I —sa br \ 
\ u{t - To) ba J 



(2.12) 



is the Jacobian of the map {a, r) — > (p, z). 

In the limit of zero velocity dispersion, the density of dark matter particles is infinite at 
the location of caustic surfaces. Thus, the location of the caustic ring surface is obtained by 
demanding that 



D2{CK, t) — —b[uT{T — To) + SO^] — 



(2.13) 



This implies 



(2.14) 



with < ;^ < 1, and hence [4] 



Pir) 



a + -{T-To){2T-To) , z(t) = ±6^-t3(to - r) , 



(2.15) 



for the location of the tricusp outline. Near the surface of a caustic ring, but away from the 
cusps, the density profile is that of a simple fold: d{a) = ^0(cr), with a > inside the 
tricusp. In the next three subsections we calculate the fold coefficient A and the principal 
curvature radii, Ri and R2, at arbitrary points on the surface of a caustic ring, other than 
the cusps. As a warm-up, we start with a special point, marked by a star in Fig. 2. In the 
fourth subsection, we obtain the density profile near the cusps. 
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1. A sample point 



As an example, we determine the fold coefficient A at (p, z) = (a, 0); see Fig. 2. Setting 
a — z — 0,a = p — a = |m(t — Tq)^ and r ~ Tq, we have 

|L>2(t)| ~26Vp^ . (2.16) 

Including the factor 2 for in and out flows, Eq. (2.11) yields: 

d(a) = — — — . (2.17 



Therefore, 

d^M 1 1 



dO^dt ah -s/p 



(2.18) 



where the subscript is to indicate that A is being evaluated at the sample point. 

To obtain a numerical estimate for Aq, we again use the self-similar infall model [12,13,11] 
with e = 0.2. For the nth ring, we have [3] 

{a„ : n = 1, 2, ...} ^ (39, 19.5, 13, 10, 8, ...)kpc ■ ( ^) f ^) f — , (2.19) 



0.27/ \ h J V220km/s^ 

where jmax is a parameter, with a specific value for each halo, which is proportional to the 
amount of angular momentum that the dark matter particles have [11]. Also, 



d?M 



dQdt 



where Vn is the velocity of the particles in the nth caustic ring, and the dimensionless 
coefficients characterize the density of the nth in and out flow. In the self-similar model 
[3] 

{/„ : n = 1, 2, ...} ~ (13, 5.5, 3.5, 2.5, 2, ...) • 10"^ (2.21) 

for e = 0.2. The /„ are like the F„ in Eq. (2.7), but they describe the nth in and out flow 
near the caustic ring whereas the describe that flow near turnaround. 
Combining Eqs. (2.18) and (2.20) we have 

A _ ^rot fn '^n 1 /r. 

^0,n — -T-^ 7 ■ {I. II) 

AttG ttn On WPn 



It was shown in ref. [4] that bn and Vn are of the same order of magnitude. Also, in ref. 
[8] , ten rises in the rotation curve of the Milky Way were interpreted as the effect of caustic 
rings. In that case, the widths p„ of caustic rings are determined from the observed widths 
of the rises. Typically one finds p„ ~ 0.1 a„. Using this and ~ J'n, Eq. (2.22) yields the 
estimates 
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{Ao,„:n = l,2,...}~(3, 4, 4, 5, 5,...^ ^ 



/0.27\^'^^ f h \ ^^'^ f Vj-ot 
'[j^J [o7l) 1,220 km/s J ■ ^^'^^^ 

At the point under consideration, {p,z) — (a, 0), the surface of the caustic ring is convex 
for all lines of sight, i.e. all tangents at that point are on the side with two extra flows. 
If the line of sight is in the z = plane, the curvature radius is a. If the line of sight is 
perpendicular to the z = plane, the curvature radius is 2-^p. Lensing by a convex caustic 
surface is discussed in subsection IV. B. 

To obtain the lensing properties of the caustic ring surface at an arbitrary point, we 
need the coefficient A and the curvature radii Ri and R2 at all locations. We derive these 
quantities in the next two subsubsections. 



2. The fold coefficient everywhere 

We choose an arbitrary point on the tricusp, i.e. on the surface of the caustic ring. Its 
parameters are (ai, ri) with ai given in terms of ri by Eq. (2.14). The physical coordinates 
(pi, Zi) are given in terms of ri by Eq. (2.15). We assume that the point is not at one of 
the three cusps. The latter are located at Ti = 0, at Ti = |to with ai > 0, and at ti = |to 
with «! < 0. 

The vanishing of D2{ti) implies the existence of a zero eigenvector of the matrix 

^(n) = ( , l"' \ ■ (2.24) 

\ U{Tx - To) hOLx } 



Let us define ^i(ti) such that 

cos {Q\) 



^'W Z ; . (2.25) 



We have 



sin (^1) = = and cos (^1) = = . (2.26) 

/(6ri)2 + {sa^Y ^(6ri)2 + (sai)'' 



We define new Cartesian coordinates (cr, r\) related to (p — pi, 2; — Z\) by a rotation of angle 
^i + f: 



^ ^ _ - sin (^1) - cos (^1) \{ P-P\\ 

We now show that a is the coordinate in the direction orthogonal to the caustic surface at 
(Pi,^i)- 

Consider small deviations about (cci, ti) in parameter space: (a, r) = (q;i + Aa, ti + At). 
Eq. (2.12) implies. 
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Ap 
Az 



At 



+ 0(Aa^ AT^ AaAr) 



(2.28) 



where T indicates transposition. The expansion of a in powers of Act and At yields: 

a = 0{Aa^, Ar^ AaAr) , (2.29) 
because the first order terms vanish: 

' sin (^i) ^ 



da 
da 



("l,Tl) 



At = -(Act At) V{ti) 



(ai,n) 



cos (6*1) 







(2.30) 



The fact that a is second order in Aa and At shows that a is the coordinate in the direction 
perpendicular to the caustic surface, and 



0- = — sin {9i)p — cos {9i)z 



(2.31) 



is the unit normal to the surface, pointing inward. 9{ti) is the angle between the p axis and 
the tangent line at {pi,Zi). See Fig. 3. 

To obtain the density profile of the caustic near the point under consideration, we need 
D2{r]i cr) to order y/o-. So we calculate cr to second order in powers of Aa and At, and D2 
and 77 to first order. We find 



and 



a--iAa AT)i '''''^^'^ -bcos{e,)\ ( Aa\ 
'^^ 2^ \-bcos{9i) -usm{9i) )[At ) ' 



D2 = -b [2saiAa + u{2ti - tq) At] 
77 = [u{ti — To) COS 9i — bai sin ^1] At — [sai cos 9i + 6ti sin ^1] Aa 



(2.32) 



(2.33) 



Eqs. (2.33) can be inverted to obtain Act and Ar as functions of D2 and rj. When the result 
is inserted into Eq. (2.32), we obtain 



'^{D2,V) 



2^(6ti)2 + (sai)2 Mri|3To - 4n\ 



2x2 



2 '2 



{bny + {sa,y 



(2.34) 



This imphes 



62 



\ l + (^-l)^ 

\ ^US 'TO 

Along the a direction {rj = 0), we have 

D2{a) = 2C{ti) b^ 

with 



us 



+ 2b'^u\3To - 4ti I ^/ — (to - Ti)ti + Tf cr 



(2.35) 



(2.36) 
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Ti , US 



1 - 



To V "^0 To 



+ 



'To' 



(2.37) 



Combining Eqs. (2.11) and (2.36), and minding the factor of two because two flows con- 
tribute, we have 



(2.38) 



with 



cos (cti) 



In terms of Ao^n = ^n(Ti = tq), for which estimates are provided in Eq. (2.23), we have 

On cosq;i(ti) 



2.39) 



pnin) C„(ti) 



2.40) 



Note that A{ti) diverges at each of the three cusps because C vanishes there. In ref. [4] the 
caustic ring parameters (a, r, b, u, s) are related to the velocity distribution of the flow at 
last turnaround. 



3. Curvature radii everywhere 

The lensing by a caustic surface depends upon the curvature radius R of the surface 
along the hne of sight. R is given by: 

1_ ^ (coso;)^ ^ {sinuf 
R R\ R2 

where Ri and R2 are the principal curvature radii of the surface and u is the angle between 
the line of sight and the direction associated with Ri. In this subsubscction, we adopt the 
convention that R is positive (negative) if, along the line of sight, the surface curves towards 
(away from) the side with two extra flows. If R is positive, the surface is called 'concave'; 
see Fig. 4. If R is negative, the surface is called 'convex'; see Fig. 6. It is a straightforward 
exercise to calculate the curvature radii of the caustic ring surface at an arbitrary point 
(q;i(ti), Ti). We find 




Ri{ti) = -2^pW- |3 - 4:^1 1 - (1 - -):^ (2.42) 

\ To To V SU TqJ 

in the direction of the cross-sectional plane of the caustic ring, and 

fl,K) = ±^ = ±^J^^T7^^(a+f(T„-r0(r.-2r0) , (2.43) 
sm6'(ri) V Ti su \ 2 J 
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in the direction perpendicular to the cross-sectional plane. In Eq. (2.43), the + sign pertains 

if < ^ < |, whereas the - sign pertains if | < < 1. Ri is always negative except at 
the three cusps where it vanishes. R2 diverges at the cusp in the z = plane. 



_ ^ 3 
"0 — 

are at angles 



For < ^ < |, there is a pair of lines of sight for which the curvature vanishes. They 



/ R2 

a; = ± arctan y - ^ (2.44) 

relative to the cross-sectional plane. Gravitational lensing by a fold of zero curvature is 
discussed in section IV. C. 



4- Density near a cusp 

In this subsubsection, we derive the dark matter density profile near a cusp. For the sake 
of convenience, we choose the cusp in the z — plane at p = a + p = po, where a — t — 0. 
Very close to the cusp, we may neglect the term of order in Eqs. (2.10). Hence 

z = bar ^ p = Po — utqt — -a^ . (2.45) 

The term of order o? cannot be neglected. We define new dimensionless quantities [4]: 

A^^,/^, T^^, X^^, Z^^, (2.46) 

To V M To V p 

where C = p- In terms of these, Eqs. (2.45) become 

Z = 2AT , (2.47) 
X = -2T - . (2.48) 

The Jacobian, Eq. (2.13), becomes: 

D2{A,T)^2hp{T - A^) . (2.49) 

Substitution of Eq. (2.47) into Eq. (2.48) yields the third order polynomial equation: 

A^ + XA + Z = 0. (2.50) 

The discriminant is: 

^-flV+f^y. (2.51) 



2/ V 3 

If 5 > 0, the cubic equation has one real root, and two complex roots which are complex 
conjugate of each other. If 5 < 0, all the roots are real and unequal. For 5 = 0, all the roots 
are real and at least two are equal. The number of real roots is the number of fiows at a given 
location. The tricusp has two flows outside and four inside. However, in the neighborhood 
of a cusp, one of the fiows of the tricusp is non-singular. It does not participate in the cusp 
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caustic. To include the root corresponding to the non-singular flow near (z, p) — (0, po), one 
must keep the term of order in Eq. (2.45). 

5 = is the equation for the caustic surface in physical space. Indeed, Eq. (2.49) implies 
that at the caustic T = A^. Therefore, Z = 2T^/'^ and X = -3T. A straightforward 
calculation shows that 

{A^ + lA^T - ST^) . (2.52) 



54 6p 

Thus = implies 5 = 0. However, the converse is not true: 5 = does not imply = 
because not all flows at the location of the caustic surface are singular. 
Eq. (2.11) for the density becomes 

1 <pm ^1 

2pohpdmtf^^\T-A''\^ ' ^ 

where the sum is over the flows, i.e. the real roots of the cubic polynomial (2.50). If 5 > 0, 
the one real root is: 

A={-^ + Vsy/' + (-| - Vsy/'. (2.54) 

It describes the one flow outside the cusp. Using Eqs. (2.48, 2.54) in Eq. (2.53) we obtain 

1 d^M 1 
d — (2 55) 

Just above or below the cusp, where X — Q and |Z| <^ 1, we have 

1 d^M 1 

2,hppodmt\Z\^l^ ■ ^ ^^^^ 

On the other hand, if we approach the cusp in the plane of the ring from the outside, where 
Z = and < X < 1, we find 

^_ 1 d'^M 1 
h po dVtdt {p — Po) 

Next, we calculate the density inside the cusp, where 6 < 0. The three real roots of the 
polynomial (2.50) are: 



A^^2\l^ cose (2.58) 



^ = 2y^cos(e + ^) (2.59) 



^3 = 2W^cos(^+y), (2.60) 



14 



where cos 3^ = (~^) ^ < ^ < |. Inserting them into Eq. (2.53) and using Eq. 
(2.48), we obtain 

~ bppodVldt V^J 1^4cos2(^) - 1 4cos2(^ + f ) - 1 ^ 1 _ 4 cos2(^ + f ) J ' ^ " 

where each term in the parentheses corresponds to one of the three flows. Adding the 
individual flow densities yields: 

2 d^^M 11 , , 

d= ■ 2.62 

bppodndt\X\{^-tein0)sm20 ^ ^ 

If we approach the cusp in the galactic plane from the inside, where Z — and < —X <^ 1, 
we have 

^_ 2 d'M 1 

b po dQdt (po — p) ' 

which is twice the result in Eq. (2.57). The gravitational lensing properties of a cusp are 
derived, for a special line of sight, in subsection IV. D. 



III. GENERAL FORMALISM 

The flrst part of this section gives a brief account of the gravitational lensing formahsm 
[14] . In the second part we show how this formalism can be streamlined in the special case 
of lensing by dark matter caustics. 

In linear approximation, the deflection angle ^ of a light ray due to a gravitational fleld 
is given by: 

^^l^dy (3.1) 

where $ is the Newtonian potential. We choose the y-axis in the direction of propagation 
of light. The deflection angle 9 is related to the angular shift ^/ — on the sky of the 
apparent direction of a source: 

m^^^-Cs) (3.2) 

where Ds and D^s ^''^^ the distances of the source to the observer and to the lens respectively. 

is the angular position of the source in the absence of the lens whereas ^/ is the angular 
position of the image with the lens present. The angles carry components in the x and 
2;-directions: 9 — {9x,9z), C — iCx,Cz), etc. Unless otherwise stated, we mean by a vector a 
quantity with components in the x- and ^-directions. 

It is convenient to introduce a 2-D potential 

or) /• 
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so that 

0- ^^V^zV^lC/) , (3.4) 
where V^^ = -DlV. -Dl is the distance of the observer to the lens. Eq. (3.2) then becomes: 

ii = is + ^^Mii) ■ (3-5) 

It gives the map ^s{^i) from the image plane to the source plane. The inverse map may 
be one to one, or one to many. In the latter case, there are multiple images and infinite 
magnification when a pair of images merges. 

Our problem is to find the image map of a point source when the matter distribution is 
given. The 2-D gravitational potential t/^ obeys the Poisson equation: 

where S(^/x, ^iz) is the column density, i.e. the integral of the volume density along the line 
of sight: 

T.{ii) = j dyd{DLii.,y,DLii,) , (3.7) 
and Ec is the critical surface density: 

It is such that a uniform sheet of density Eg focuses radiation from the source to the observer. 
Eq. (3.6) is solved by: 

V'(6) = ^/ d2e;S(g)ln|6-g| , (3.9) 

and hence 

AC = ii - is = = / d% . (3.10) 

The image structure, distortion and magnification are given by the Jacobian matrix of the 
map CsiCi) from image to source: 

where ipij = Q^-^p. Because gravitational lensing does not change surface brightness, the 
magnification A4 is the ratio of image area to source area. Therefore 

^^^^ ■ 
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To first order, for ijjij <ti 1, 

M = l + Vlip = l + 2^ . (3.13) 

To get the largest lensing effects, we wish to minimize given in Eq. (3.8). For fixed Ds, 
the minimum occurs when the lens is situated half way between the source and the observer. 

Also Ds should be as large as possible. For our estimates, we will assume that the source is 
at cosmological distances, e.g. 2Dl = 2Dls = Ds = IGpc, in which case Sc = 1.39 g/cm^. 

For a general mass distribution, the gravitational lensing effects are obtained by first 
calculating the column density, Eq. (3.7), and then the image shift, Eq. (3.10). This 
procedure can be simplified when the lens is a dark matter caustic. We are interested in 
lines of sight which are tangent to a caustic surface, because the column density S has the 
highest contrast there. We assume that, in the neighborhood of the tangent point, the flow 
is independent of y except for a shift x{y) of the caustic surface with y. The density can 
then be written as: 

d{x,y,z) = d{x - x{y),z - z{y)) , (3.14) 

where d{x, z) is the density of the 2-D flow in the plane orthogonal to the line of sight. The 
flow in the y-direction is irrelevant because lensing depends only on the column density E. 

In the limit of zero velocity dispersion, a 2-D flow is specifled by giving the spatial 
coordinates x{a,(3,t) of the particle labeled {a, (3) at time t, for all {a,(3,t). The labels a 
and /3 are chosen arbitrarily. At a given time t, x{a, /3,t) — x has a discrete number of 
solutions {a, f3)j{x,t) labeled by j = 1, ...,n{x,t). n is the number of distinct flows at {x,t). 
Here, t is the time at which the light ray passes by the caustic. Henceforth we will not show 
the time dependence explicitly. The particle density in physical space is: 



"^^^ d^A 1 



(3.15) 



where A is the mass per unit length in the direction of the line of sight and is the A 



^ l9{a,/3) 

dadfz 

density in parameter space. Inserting Eqs. (3.7) and (3.14) into Eq. (3.10) we obtain: 

— * — * 

When d is replaced by the caustic density (3.15), Eq.(3.16) becomes: 



1 f ^ x-x' ^ 1 d^A 



(3.17) 



{a,(3)={a,l3)Ax'-x{y)) 



Changing variables from {x', z') to (a, P) and assuming that the density in parameter space 
varies only slowly over the region of integration, we rewrite Eq. (3.17) as: 

V„«K, = x/BJ = — — . (3.18) 
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Further simplification is achieved by defining the complex integral: 

I(x,z) = [ dad(3 -. — — — \-. -. — — - , (3.19) 

in terms of which the shifts are given by: 

X ^ dip 1 (fA ^ f , ^, 

= tt:;— = ^ ^ - — — Ke ay lix — xiy)) , 

Eqs. (3.19) and (3.20) are useful when the caustic has contrast in the two dimensions 
transverse to the hne of sight, e.g. near a cusp. 

However, in many apphcations, the caustic has contrast in only one of the dimensions 
transverse to the hne of sight. Choosing x to be the trivial direction, Eq. (3.6) is reduced 
to: 

^(6) -^n^- D,Ci) , (3.21) 
and the column density is given by: 

mi) ^ j dy d{z - z{y)) . (3.22) 
The fiow is now effectively one dimensional. Its physical space density is given by 



"'^^ dh 1 

.7=1 I da I 



(3.23) 

a=aj{z) 



where A is the mass surface density in the two trivial directions [x and y) and ^ is the A 
density in parameter space. The 1-D Green's function is G = |(|^| + a^) + h, where a and h 
are arbitrary constants. The shift is: 

^(6) = ^Jdi'i m'l) (Sign(6 - O + a) . (3.24) 

The constant a causes an overall shift of the image, which does not concern us. We choose 
a — —1. Repeating the steps done earher for the 2-D case, Eq. (3.24) is re-expressed as: 

2 dA r r 

AC = J dyj da Q{-z + z{a) + z{y)) . (3.25) 

In the next section Eq. (3.25) is used to calculate the shifts due to simple folds of dark 
matter fiows, and Eq. (3.20) for the shifts due to a cusp. 
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IV. APPLICATIONS 



A. Lensing by a concave fold 

We consider a simple fold caustic which has curvature radius R along the line of sight. 
We assume that the surface is concave, i.e. it is curved in the direction of the two extra 
flows. See Fig. 4. The outer caustics of dark matter halos are examples of concave caustic 
surfaces. The convex case is discussed in the next subsection. 

In the neighborhood of the point P of closest approach of the hue of sight with the 
caustic surface, we choose coordinates such that z is perpendicular to the surface whereas 
X and the direction y of the line of sight are parallel. P has coordinates x — y = z = In 
the 1/ = plane, the flow is given by z{a) = —-^ha"^, where h is a positive constant. Using 
Eq. (3.23) we find the density in the y = plane: 

dia) = A^ , (4.1) 



where a — —z, and 



/2 dA 



2 

For 1/ 7^ 0, the density is still given by Eq. (4.1) but with a — z{y) — z and z{y) — — f^- We 
calculate the shift using Eq. (3.25): 

2 !^ ...i^Ql_,h._y'\ 



AC^-——— / dy / . 



-'c 

AndA IR 



2 2R 



Hence 

A^ = ^i-^s = V^iO{-^i) , (4.4) 

with 



V = ■ (4.5) 



One can also obtain this result by calculating the column density: 



(4.6) 



= [dy d{y, z = Ci)= [dy / ' 2 = ^AV2Re{-ii) , 

and solving Eq. (3.6): 

ii-is=^^^^J di: S(6) = r,i: e(-6) • (4.7) 
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Eqs. (4.6, 4.7) were first obtained by Hogan [6]. The agreement with Eq. (4.4) vahdates 
the formahsm derived in section III. Fig. 5 plots the source position ^5 versus the image 
position 

The magnification is: 

M = ^^i + vei-Ci) + oiv') ■ (4.8) 

When the hne of sight of a moving source crosses the surface of a simple concave fold, 
the component of its apparent velocity perpendicular to the fold changes abruptly. Also, a 
discontinuity occurs in the magnification of the image. Both effects are of order rj. 
We estimate rj for the outer caustics of galactic halos using Eqs. (2.5-2.8): 



V2 



K = : n = 1, 2, ...} ~ (7, 6, 6, 6, 6, ...) • IQ-^ • 

fDLDLs]( Vro. \(h_\ .^g. 

\Ds Gpc J V220km/s J \0.7 J ' ^ ' ^ 

A magification of order 10~^ seems hard to observe. However, the images of extended sources 
may be modified in recognizable ways. In particular, straight jets would be seen with an 
abrupt bend where their line of sight crosses a fold. Indeed the image is stretched by the 
factor 1 + 77 in the direction perpendicular to the caustic, on the side with the two extra 
flows. If the jet makes angle a with the normal, it appears bent by an angle 5 = ^-q sin {2a) . 
Searching the sky for bends in extended sources may be a realistic approach to detecting 
caustic structure in galactic halos [6]. 



B. Lensing by a convex fold 

By deflnition, a convex fold is curved in the direction opposite to the side with two extra 
flows. See Fig. 6. Using the conventions of the previous subsection, we write the equation 

2 

for the displacement of the surface along the line of sight as z{y) = |^ , and that for the 
flow as z{a) = —^ha"^ , for smaU z, with R,h>0. Eqs. (4.1) and (4.2) still apply, with 
(7 = y^/2R - ^ . Eq. (3.25) yields the shift: 



Ae = -^;^ e(-e,)/ dy + Q{ij) dy\ 
T.cDl da y Jo J^pRiHii J 



(4.10) 



We introduce a cut-off L for the integral over large y. L can be thought of as the length 
scale beyond which our description of the flow is invalid. The above equation becomes: 



AC 



4 dK 1 
Ec-Dl da ^/hR 



LJL^ - 2RDl^i - 2RDl^i In 



L + VL^ - 2RDl^i 
pRDL\b\ 



■ (4.11) 



By expanding Eq. (4.11) in powers of and using Eq. (4.2), we obtain the ^/-dependent 
shift: 
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7] 



In 



TT 



where rj' is [like rj in Eq. (4.5)] given by 



2nAx/2R 



rj = 



6 



(4.12) 



(4.13) 



The magnification is : 



M = 



1 + ^ In 



TT 



2L2 



(4.14) 



The cut-ofT L has an effect on both the magnification and the elongation of the image in 
the direction normal to the caustic surface, but that effect is ^/ independent. L has a global 
effect on the image, as opposed to an effect locafized near ^/ = 0. 

Fig. 7 plots ^5 versus It shows that a convex fold can cause a triple image of a point 
source. In particular, when the source is exactly behind the caustic (^5 = 0), the images are 
at — — 0, and +^c, with 



. ^ TT 



(4.15) 



Sufficiently far from the caustic, a point source has a single image, say at > 0. When the 
line of sight approaches the caustic surface tangent point, two new images appear on top of 
each other at = $,13 — — ^c/e. At that moment, the magnification at $,12 is infinite, and 

^5 = Vj^f-je-n. As the source crosses the caustic, moves towards ^/i and finally merges 
with it. When ^/i = (^/2 = +Cc/e, the magnification diverges again. After that, only the 
image at ^/a remains. 

Let us apply the above results to the line of sight in the plane of a caustic ring at the 
sample point, (z = 0,p = a), discussed in subsection II.B.l. Setting R — an and using Eq. 
(2.22) we obtain for the nth ring 



/ 27r Ao,n"\/2a^ 



^rot 



Using Eq. (2.23) to estimate /lo,n, we find 

{77; 1,2,...} -(5, 4, 4, 4, 4,...) -10 



(4.16) 



2 Dls 



'0.27 

I Jmax . 



h 



■^rot 



^0.7y V220km/s^ 



(4.17) 



For such small values of rj' the angular distance between the triple images is exponentially 
small and unresolvable with present and foreseeable instruments. 

The image of an extended object is stretched in the direction perpendicular to the caustic 
by the relative amount 
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. (4.18, 



TT 

The image is stretched for ^7 < and compressed for ^/ > where 

2L2 



RD, 



(4.19) 



In the case of the sample point {z = 0, p = a) with hne of sight in the z = plane, the cut-off 
length L, i.e. the distance in the |/-direction over which our desciption of the flow is valid, 
is of order y/2ap. In that case ~ ^p/Dl. Since p/Dl is the transverse angular size of 
the caustic ring, our description certainly fails for ^7 > p/D^- So, over the region where our 
calculation is valid, the image is magnified. The effects are generically of order one percent, 
increasing to several percent when ^7 ~ 10"^^^, for caustic rings at cosmological distances. 



C. Lensing by a fold with zero curvature 



We saw in section II. B. 3 that the surface of a caustic ring has tangent hues along which 
the curvature vanishes. One may speculate that the lensing effects of a caustic surface are 
strongest when the line of sight is tangent to the surface in a direction of zero curvature, 
because the line of sight stays close to the caustic over greater depths in that case. If the 
line of sight of an observer looking at the outside profile of a caustic ring is, at some point 
on the profile, in a direction of zero curvature, then the equation for the intersection of the 
caustic surface with the plane containing the outward normal to the surface (i) and the line 
of sight {y) is ziij) = where U is positive and has dimensions of (length)^. In such 
a case, the cubic term in the Taylor expansion of z{y) is absent because the line of sight 
remains everywhere outside the caustic ring tube. We did not calculate U for caustic rings 
but expect U ~ (kpc)"^ in order of magnitude. The flow is given by z{a) = —\ha^ with 
/i > 0, as before, and Eq. (4.2) holds. 

Using the above expressions for z{y) and z{a) in Eq. (3.25), we find the shift 



AC 



2 dk 



da Q \ —z ha^ — 

\ 2 4C/ 



Sc-Dl da 
J^cDida ^/2M7 ^ ^ Jo 

--e(-6)(-eoel)''' , 



where 



(4.20) 



(4.21) 



The magnification is: 



M = 



1/4 



(4.22) 
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Fig. 8 shows vcrsTis ^j. Triple images occur when |,^/| ^ ^o- Unfortunately, for the zero 
curvature tangents of caustic rings envisaged above, ,^0 is very small. For 2Dl = 2D is = 
Ds = Gpc, A = 3 - 10-'^ 1 , and f/ = (kpc)^ one finds 6 = 4- 10"^^. Hence the 

cm-^kpc^ 

triple images cannot be resolved. Also, even at angular distances as small as ^7 ~ 10~^, the 
magnification and image distortion is only of order 1%. 



D. Lensing by a cusp 

In this subsection, we investigate a line of sight parallel to the plane {z = 0) of a caustic 
ring, and passing near the cusp at p = po- See Fig. 9. We use the 2-D lensing equations 
derived in section 111. The shifts are given by Eq. (3.20) in terms of the complex integral / 
ofEq. (3.19). Using Eq. (2.45), we have 

7 = [da [ dr 1 — 2 — 5 7 ■ (4-23) 

J J p — a — -^utq -\- utqt -\- -^sa'^ -\- lyz — bar) 

Because we are close to the cusp, the term of order is neglected in the denominator of 
the integrand. In terms of the parameters defined in Eqs. (2.46), 

2 /■ /-^ dA 
^^W^J Loo A^-2i^ + 2T + X + t^ ■ ^^-^^^ 

The integration over A yields 

J _ 2Tr r (IT 



where C = Re + 2T + X + > . Near the cusp (X,Z « 1) the terms of order 

can be neglected since they are unimportant in the denominator of Eq. (4.24). Eq. (4.25) 
becomes then: 

IiX,Z)^^f , (4.26) 
WQJ l2T + X + i^ 

where the integration domain is defined by the inequality: 

^4 _ ^^2^3 + XT^) - < . (4.27) 

Let us call Tj [i = 1, ...,4) the roots of the polynomial on the left hand side of (4.27). Near 
the cusp, three of the roots are near zero and one of the roots is close to 2^. Let us call 
the latter T4. In order to find the roots near T = we neglect the quartic term and solve 
the cubic equation: 

X Z^ 

+ — + — = . (4.28) 
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This equation is also obtained by eliminating A from Eqs. (2.47) and (2.48). Hence the 
solutions of Eq. (4.28) are the values of T for the flows near the cusp. As was discussed in 
section II. B. 4 , the number of flows is determined by 6, given in Eq. (2.51). Outside the 
cusp, where 5 > 0, there is a single root Ti. Inside the cusp where S < 0, there are three 
roots Ti, T2 and T3. Fig. 10 shows the quartic polynomial, Eq. (4.27), for various values of 
X and Z, and C = 1- 

For a line of sight just outside the cusp, the quartic polynomial is negative from Ti to 
T4 with: 

Eq. (4.26) becomes in that case: 



hr^ + X + i^- J2Ti + X + i^\ . (4.30) 



bVCi^ ^ VC V ^ vc 

Let us consider the hne of sight defined by Fig. 9. Eq. (3.20) gives the shift in the direction 
perpendicular to the plane of the cusp as 

2 

where X(y) = —2fm> shift of the cusp as a function of depth. The integral in Eq. 

(4.31) with the integrand given by Eqs. (4.29) and (4.30) can be done numerically. Here we 
only give a rough estimate, to determine the order of magnitude and qualitative properties 
of the lensing effects. 

For X — and Z « 1, we have 



-j=Sign(Z)|Z|V3_^|Z|i 



+ 0{Z) . (4.32) 



Because the scaling law X^ ~ Z^ holds close to the cusp, we expect Eq. (4.32) to be valid as 
long as the shift in the x-direction |AX| ^ \Z\^. Hence Eq. (4.32) provides a good estimate 
of I{X — X{y), Z) over a depth 2 Ay with Ay ~ \/2p^ |Z| 3 . Therefore 



Ae. S,g„(Z) . (4.33) 



HcDhb V C dadT 
Since we are near a = 0, 

d'^A 1 d^M 1 d'^M 



(4.34) 



dadr 27rpo dadr po dVldt 
Using this and Eq. (2.18), we obtain 

Ae.~r7"|6.|^Sign(e/.) , (4.35) 
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where 



2V2a 



The contribution to the magnification from distortion in the ^-direction is therefore: 



(4.36) 



(4.37) 



Fig. 11 plots ^sz versus ^/^. For ^/^ ^ ^1 ~ i] , a point source has triple images. For the 
nth caustic ring, 77" can be estimated using Eqs. (2.19, 2.23), and assuming p ~ 0.1a, C ~ 1- 

-4 Dl Dls 



{7?::n = l,2,...}~(1.5, 1, 0.7, 0.6, 0.5 




(4.38) 



At cosmological distances, the typical angular separation between the images is of order 
10~^^. Unfortunately, here again, the images are too close to be resolved. However, for an 
angular distance ^7 ~ 10~^, M.z is of order 10%, which may be observable. The lensing 
properties of a cusp for other hues of sight can be calculated numerically using Eqs. (3.20) 
and (4.23). The latter equation is replaced by Eqs. (4.30) and (4.29), and 74 = 2C, if the 
line of sight is everywhere outside the region with two extra flows. 



V. CONCLUSIONS 

The dark matter caustics present in the halos of isolated galaxies were discussed in 
Section II. They are of two types, outer and inner. The outer caustics are topological spheres 
surrounding the galaxy and labeled by n = 1,2,3... Each is a simple fold (^2) catastrophe 
located near where a dark matter outflow reaches its furthest distance from the galactic 
center before falling back in. The inner caustics are rings in or near the galactic plane. 
They are also labeled by n = 1, 2, 3... A caustic ring is a closed tube whose cross-section is 
an elliptic umhilic {D_i) catastrophe. The cross-section, shown in Fig. 2, is characterized 
by three cusps. A caustic ring is located near where the particles with the most angular 
momentum in a given inflow are at their distance of closest approach to the galactic center 
before falling back out. In section II, we gave estimates of the radii and the fold coefficients 
of the outer caustics, based upon the self-similar infall model. For the inner caustic rings, 
the curvature radii and fold coefficients depend on position on the surface of the ring caustic. 
We derived formulae for these quantities as a function of position (parametrized by Ti) and 
of the five caustic ring parameters (a, 6, m, tq, s). A great fraction of the surface of a caustic 
ring (0 < 1^1 < |) is saddle-shaped and therefore has tangent directions along which the 
curvature of the surface vanishes. Finally, in section II, we derived the mass density profile 
near a cusp caustic. 

Gravitational lensing was discussed in section III. In general, the shift in image position 
is the gradient of a potential whose 2-dim Laplacian is the column density. For an arbitrary 
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mass distribution, the procedure for calculating the shift involves two steps. First the 
matter density is integrated along the line of sight to obtain the matter density. Second, the 
potential is obtained by convoluting the column density with the 2-dim. Green's function. 
In the case of lensing by dark matter caustics, this procedure can be simphfied: the shift is 
expressed directly as an integral over the parameter space of the dark matter flow forming 
the caustic. The relevant result is given in Eqs. (3.19) and (3.20) if the caustic has contrast 
in the two dimensions transverse to the line of sight, and in Eq. (3.25) if the caustic has 
contrast in only one of the dimensions transverse to the line of sight. 

In section IV, we applied our formalism to the gravitational lensing by dark matter 
caustics in four specific cases. In the first three, the hue of sight is tangent to a surface where 
a simple fold catastrophe is located. The three cases are distinguished by the curvature of 
the caustic surface at the tangent point in the direction of the line of sight: 1) the surface 
curves toward the side with two extra flows, 2) the surface curves away from the side with 
two extra flows, and 3) the surface has zero curvature. In the fourth case (4) studied, the 
hne of sight is at a specific location close to a cusp catastrophe and parallel to the plane 
of the cusp. We found that each case has characteristic lensing signatures. In three of the 
cases (2, 3 and 4) there are multiple images and infinite magnification of point sources when 
their images merge. In case 1, there are no multiple images. Unfortunately, the effects are 
small even for dark matter caustic lenses at cosmological distances. The multiple images 
of a point source cannot be resolved with present instruments. Typical magnifications and 
image distortions are of order one % to a few %. 
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FIG. 1. A snapshot of the phase-space distribution of CDM particles in a galactic halo. For 
simplicity, spherical symmetry is assumed, r is galactocentric distance and v is radial velocity. 
The solid line indicates the location of the particles. The dotted line corresponds to the observer 
position. Each intersection of the solid and dotted lines corresponds to a CDM flow at the observer's 
location. "Turnaround" refers to the moments in a particle's history when it has zero radial velocity 
with respect to the galactic center. A "caustic" appears wherever the phase-space line folds back. 
Particles pile up and hence the density diverges there. 
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FIG. 2. Cross-section of a caustic ring in the case of axial and reflection symmetry, p and q 
are the transverse dimensions in the p and z directions respectively, a is the ring radius. The 
star indicates the sample location (p ~ a, z = 0) discussed in section II.B.l. For clarity, we took 
p,q ^ a. For actual caustic rings, p,q « a. 




FIG. 3. An arbitrary point on the tricusp is labeled by ri with ai = ai(Ti). Its physical 
coordinates are {pi,zi). A new Cartesian coordinate system (cj, r/) is defined there such that a is 
perpendicular to the caustic surface. It is rotated relative to the (p, z) cooordinates by an angle 
e{n) + f . 
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FIG. 4. Lensing by a concave fold. The arc is the intersection of the caustic surface with the 
plane containing the normal (i) to the surface and the line of sight (y) . The shaded area indicates 
the side with the two extra flows. 
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FIG. 5. Source position ^5 as a function of image position ^/ for lensing by a concave fold, ij is 
given by Eq. (4.5) in terms of the fold coefficient and the curvature radius of the caustic surface. 
Estimates of rj for the outer caustics of galactic halos are given in Eq. (4.9). 
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FIG. 7. Source position ^5 as a function of image position ^7 for lensing by a convex fold, is 
given in Eq. (4.15). 
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FIG. 8. Source position ,^5 as a function of image position ^7 for lensing by a fold with zero 
curvature. ^0 is given in Eq. (4.21). 
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FIG. 9. Lcnsing by the caustic ring cusp at (z, p) = (0, po), with hue of sight parallel to the 
= plane. We define x = p — pQ. a) Side view in the direction of the line of sight. The latter is 
presented by the dot near x = z = 0. b) Top view. The curve is the location of the cusp in the 
= plane. 
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